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X* Introduction 

Th« task of providing an adequate coverage for a space mission of 
the extent of tba Apollo project lsamdiately lays bare the striking defi- 
claaoaa of tha existing tracking and orbit prediction programs in use to- 
day* Cm vary nature of the Apollo mission, with its near earth circular 
parking oTblt, its highly elliptic tranalunar trajectory, its lunar capture, 
xeadesrous and earth return logs, and, finally, its critical re-entry phase, 
all require individual attention. Perhaps the most stringent requirement 
la the entire mission la tha need for precision predictions of the orbit 
over its entire 7-lV day period. Shis requirement alone is sufficient to 
eliminate most programs presently using the conventional Coveil's integration 
proced u re. In addition, tha affect of small changes over large prediction 
time area calls for a choice of differential^ correction parameters which do 
not lose information due to their secular character. Since the information 
content of tha observation over different portions of the trajectories vlll 
vary greatly vith tha vehicle position, it also becomes necessary to absorb 
and utilise information coming from inhomogeneous observations of many 
different types, such as range, range rate, azimuth and elevation, optical 
and other instnaMnt readings mads in mid-course. Furthermore, the need 
for real time decision processes requires the ability for rapidly carrying 
out Iteration procedures vbich converge rapidly. In particular, the con- 
ventional least square techniques now in cannon use may be taxed beyond 
their limits of linear assumptions. Thus, the requirements are to develop 
a general purpose real time tracking and orbit prediction program with the 
following features: 
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1 . A rapid, basic precision orbit prediction program with high 
numerical accuracy not subject to numerical round, off error over long 
time periods. 

2. She capability of handling circular, elliptical and hyperbolic 
orbits with equal flexibility under the action of perturbation forces. 
Special treatment will be necessary to handle the extremely rare cases 
of exactly rectilinear and parabolic cases. 

3 . She ability to absorb and utilize inhomogeneous observations 
from many different sources with equal facility. In particular, this re* 
quires forming matrices of partial derivatives of the observations with 
respect to the variational parameters rapidly. 

k, A theoretically sound statistical process for obtaining in- 
formation from observations rapidly vithout requiring large stretches of 
observations over long time arcs which unduly tax the linear assumption 
of the leset square techniques. 

5. A program combining these features is a necessary prerequisite 
for adequate ground coverage of the Apollo mission and moreover provides 
a ready and existing tracking program for almost any existing space mission. 
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fetation 


ZZ. 

t Vehicle position vector 

r Magnitude of position rector, r « (x 8 + y* + s 2 )^ 

F Pert* Men tion aoosleratica vector 

U Mess per— ter 

v 0 Velocity of vehicle 

| P er t ur bation (lleplanmasnt vector 

t, t ^ t r Present tine. Initial tine, tine of rectification, respectively 

a fled enjor axis 

a Souentrlelty of orbit 


1 

Q 

• 

Q 


Ihcll nation of ocbltal plane 
Right asoenslon of ascending node 
Argewnt of perigee 



felt vectors in the orbit plane directed, respectively, toward 
perigee and 90° frea perigee In the direction of aotlcn, aa 


need in 8ection ZZZ. 

Pjj, Py> p s Oonponente of vector P 

Qjg, <ly, Ccnpone n ts of vector Q 


B Bcoen trl c anonaly as used in Section HI 

B Elevation angle aa uaed In Section VZH 

n Mean notion a a uaed in Section HI 

d R • R 

0 0 0 

e • 

f, g, f, g Coordinate funetlona In the Kepler problem 

• ■I *1 Incremental eccentric anomaly 
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ufe Angular rate of the earth 1 a rotation 

1 , m, n Mini track direction coaines, Eq. (69) 

X Ml i y M, » * M 1 Topocentrlc coordinates of the vehicle In the local horizon, 

local vertical system vith x * 1 * positive aputh, y* * ' positive 
east, and z' ' ' positive upward along the local vertical. 
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0(t, t) Point trMM faaMtlca aatrix of the partial derivatives of t}>e 
state variables vith raapaot to tha orbit parameters 
0(^ t.y State transition aatrix of tte octet parameters 
Kt) Matrix of partial derivatives of tea Observations vlth raapaot 

to tea octet para— tun 

Q(t) Oovarla&oe aatrlx of tea ortlt paraaiatara 


Milt) Deviations la tea octet parameters 




Uaaar aatlaatnr for tea orbit paraaatara 
tkilt vaotora defined tgr Figure 3 
Ap p ar e nt range aa defined la Saetlon X 
Martial dlataaoaa to observer and telde, respectively, as 
naad la Section X 

Sartas of rafraotlon and average index of refraction as used In 


8aetioa X 


a Macon speed of light 

▼ Signal velocity In medium aa uaad in Section X 

a Path angle of signal, see Pig* 7 

Xq Qround station eight angle. Pig. 7 

f- Frequency of emitted signal as uaad In 8ectlon X 

ff Dlffaranea in frequency oatvaan emitted and received signals 

h Q , h^ Altitude of obaerrer an! of vehicle, respectively 
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Subscript* 


rabid* 

tb 

1 attracting body 

r*f«r*ao* body as us«d In 8*otlon ZZZ 

KapLar, i.#. # quantity obtained from solution of tiro-feody problsn 

value at Initial tine 

Talus at tlae of rectification 


Superscript 


XXX* Precision Orbit Prediction 


A thorough analysis ccagaring the various special perturbation 
■sthods In uaa today tor orbit prediction has been carried out ^ 

BtfWrenoe 1* Tbs general conclusions drawn in this study say be stated 
•lagly as follows t 

!• Ocrell*s aethod, w hile s ingle to prograa, consumes larger 
n s chine ccagutlng tines and Is subject to an unavoidable accuaulaticn 
of the round off error* 

8* 9ra e l l a d n a t ion of truncation and round off error cob be 
acccngllahed through either the variation of paraaetera or the »»»%• 
aethod. 

3* ®*e preference for the Ebeke aethod over the variation of para- 
coeea fSraa the simplicity of tbe equations and a great reduction 
1* confuting tine* Thus, it la possible to generate a precision program 
usi ng the &V;he aethod w hi c h will produce a solution of the equations 
of notion as precisely as required In shorter computing tine than any 
other available nethod. 

As la well knom, the Backs aethod solves tbe best local two body 
problen and Integrates the deviation froa this nocdnal trajectory. Since 
the round off error occurs only In the Integrated position, it is possible 
to ellalnpte this difficulty froa the answer by periodic orbit rectification. 



Addition! difficult!., of the conventional &>cke method, such a. numerical 
inaccuracies for circular orbit., etc., have been eliminated by using a 

•mutt® of the Kepler problem in ten. of the Initial portion and velocit 
vector*. - 

Sanations of Motion 

In a Bevtooian sys tea, the equation, of action of a particle in the 
gravitational field of n attracting bodies and subject to other perturbing 

acceleration, such as thrust, dreg, oblaten..., radiation pressure, etc. 
are given by 





*“*• tVittiau u. put into obwmUs fo™ by rof.rrlog tb« to 

body c(. 4. lb. .quatlau of .otlcn of tb. «f.«„c. boOy 



1-1 

l#m 



( 2 ) 


Subtraction of Bq. (2) Item B> (1) results In ttj equation of motion 
in the vehicle with respect to the xei-rence body „• 


H vc “ * (u T + * c ) 



r 3 

TC 


- L 

1-1 

i*m 




( 3 ) 
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Si 


B* jhthgd of Integration 

® ^3) ** latsgrmtsd directly by sew numerical scheme , there 

nnltl| after e umber of atejHby^step integrations , an accumulation ol‘ 
•wear KMch leads to inaccurate results. To avoid this loss in precisicn, it 
la c o n ve ni ent to write 1 *. ( 3 ) in the form 


Kc-K ♦« 

velocity and diaplamaant vectors can be written 


(4a) 


* • A 


(4b) 

(4c) 


»«f»»woa body is chosen so as to minimize the perturbations, i.e., 
the ona in whose sphere of influence the vehicle travels. 

IB this met hod S is taken as 


\ ■ * K ♦ u c> h 

*k 


(5) 


? ■ * (u. 


[X v 

D 

- > “1 

i-i 


. R <* 


X vi 

ci 


L 



*£* j <« 


Bq. (3) !• solved as 


described below; Eq. ( 6 ) is Integrated numerically. 
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C . Solution of the Kepler Problem 


1. Classical Solution 


The two- body orbit which results frap the solution of Eq. (5) with 
the initial conditions: 

\ (to) - l^ c (to) * fKto) 


\ (to) * ^(to) ■ R (to) 


( 7 ) 


nay be represented by many different sets of "elements", such as the semi* 
■ajar axis, a, the eccentricity, e, the inclination, 1, the right 
ascension at ascending node, ft, and the argument of perigee, and a time 
variable such as the time of perigee passage or the time of nodal crossing, 
to locate the vehicle in its orbit* A difficulty arises with this particular 
set if orbits of zero or small inclination are to be used, in which case 
the right ascension of the ascending code is not defined or poorly defined. 
Further, for orbits of zero or small eccentricity the argument of perigee 
is ill defined. Hie f.'.rst of these difficulties is removed by replacing 1, 
x and »• by two unit vectors P and Q in the plane of the motion, P directed 
toward perigee and Q at right angles to it. 

‘Ifce expressions connecting i, <*>, and ft with P and Q are determined 
by a tranr formation of the coordinate system of the orbital plane to a 
geocentric coordinate system by means of three successive rotations ( Eulerlan 
angles). Ihis results in: 


P x ■ cos 0 cos ® - sin 0 sis » cos 1 

Py ■ sin 0 cos <s ♦ cos 0 sis « cos 1 

p - sis « sis 1 

(e; 

- - cos 0 sis ® - sis 0 cos « cos 1 

fy " - sis A sis 01 + cos 0 cos ® cos 1 

Qg « COS « sis 1 

I 

In tsxas of tbs P sad Q tbs solutloo of the two-body problem Is giv»n 
s fusctlon of tbs scosntrlc sa cosily by 

\ - \ P ♦ P k Q 


\ * \ * ♦ 7 k Q 
Z k i t f sre gives by: 

m 

Xjj ■ o (cos B - e) 
y k • s^ 1 • e 2 sis E 
r k - * (l - e cos B) 


i . . JuS 
* r. 


sis B 


1 - e ) 


cos B 


Elliptic 


. n 


for the ellipse sad 


« a (coeb E' - e) 

3 ^ » - 1 sinh E 

p. ■ * (1 • • cosh fi) 

£ - - slab E 

r k 

y - ^ - C cosh E 

k r k 

for the hyperbola. 


Hyperbolic 


(9) 


The eccentric anomaly E is obtained as a function of tine 
E - e sin E - (Eo - e sin Eo) - u (t - to) 
e sink E - E - (e sinh Eq - Eo) ■ n (t - to) 


fron 

Elliptic 

(ID) 

Ejrperboilc 


The elesents a, e, P, Q are given in terms of the initial conditions by 

*o »(Bo . Ro)^ 

do m Bo • R o 
2 


(ID 


▼o" ■ Ro • Ro 


2 


-1 


-(l-SC) 

*0 u 


12 - 


< 



J 

••ji 

• ila Mo ■ 


• cos Mo - u 


- i-i-fa 


Elliptic (U) 


"*1ra i*/« 

• slab to * .. 

Vnli 

s cosh ^ ■ 1 - 

.. [<!-?)*♦ &] * 


- 1 


Hyperbolic 


s' 


P - 


*>- I s 

r o v u 


sin Eo Ho 


7 1 . «* Q - IgJQ Ho (co * *> - e) Ro 


Elliptic 


P - 


cosh fc 
ro 


Ho ■If slab So fio 


t 


f5TT Q ■ Ro ( co#il - e) $o 


Hyperbolic 


While this foraulatlon of the solution of the tvo-body problem Is 

attract ire frca a computational point of rierv, the scheae will hare 

» 

difficulties with circular and near circular orbits (since P and Q 
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are ill defined) , a somewhat paradoxial situation. However, these dif- 
ficulties be completely removed by reformulating the solution to be 
independent of the vectors P and Q and the eccentric ananaly £• 

2. Development of HASA Two Body Equations 

It is noticed that Bqs. ( 9 ) give end as linear combinations 
of P and Q and Et^s. (11) give P and Q In turn as linear combinations of 
Bo and Ac. It Is dear, therefore, that and R^ may be expressed as 
linear combinations of Ifc and Bo* 


\ - f Bo ♦ g Ro 

( 12 ) 

• • • • 

Rjt - f Bo ♦ g Ro 



f * ~ £(coeh • - 1)J * 


C ifcslnfcO-^S (cosh 0 - D - - Jil£ 


7 - - (cosh 0 - D ♦ 7® cosh 0 - - sinh 0 

a • 


♦ t-to 


f ■ slab 0 

r*a 


C - ? (cosh 0 • 1) ♦ 1 

T 


( 15 ) 


hyperbolic 


A similar technique applied to Kepler' ■ Sgs. (10) results in the folbwlng 

J 

equation, vhlch fhmlshes 0 as a function of the time* 


n(t*to) ■ ® - sin 0 ♦ sin 0 + (1 - cos 0) 

n(t-to) ■ siahO - 0* — ® slnh 0 ♦ (cosh 0-1) 

* J=ua 


Elliptic 

(14) 

hyperbolic 


IT the function f * , fa, fa, f* are defined as 




Ellipse 

hyperbola 

*x 

(«) 

» 0 - sin 0 

■ slnh 0-0 

fa 

(0) 

» 1 - cos 0 

- cosh 6-1 

fa 

(•) 

■ sin 0 

m slnh 0 

fe 

( 0 ) 

- cos 0 

- cosh 0 
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B»e solution of the tvo-body problem for both elliptic and hyperbolic orbits 
is given by: 

rJsL 


f ■ ro f a ♦ 1 

g ■ • ” fi ♦ (t - to) 


i5i - * * ft * -fist* 


f m m r&. i Jfti - £3 

i |a| r 0 r * 3 


( 16 ) 


. » 1*1 

g - r fa + 1 


n(t-to) “ fi fa * . 

1 ' /»M 


D. numerical Procedures 

It will be noted that the only distinction between elliptic and 
hyperbolic paths Is completely contained In the definitions of f 4 , fa, 
fa# f« (Bis. 1 ^). If these functions are considered as power series 
expansions in d 2 as given in Eqs . (17) , then the switch from elliptic 
to hyperbolic cases is achieved by replacing -8 2 by 6 s , permitting 
compact and efficient programming. 

Considering the first lfc terms of a power series expansion for 
sin 9 , cos 6, sinh 6 and cosh 9 , the power series expansion for ft 
and f 2 any be expressed as: 
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f x - {((O^g ♦ 1) 25?24 + 1) ••• ♦ 1) 375} 

f*«i-co«e-. C(((^Jg$ ♦ D ^753 + D ••• ♦ 1)3^} 

f a . mu e - • . 0 {(((^ ..+ 3)f 2 } 

fa ■ eoah $ * 1 • (((^ g^ * 1) ••• ♦ ') 


Elliptic 


( 17 ) 


Hyperbolic 


She UN of Cbebymher polynomials allows the calculation of t x end 
fa to the ease accuracy vith fever terns. They are currently being 



Zn order to inure a loss of accuracy, the aethod of 

confutation of f x , f a , fa, f 4 will depend upon the oegnitude of 0. 
She following testa are nade to deternlss the aethod of computation 
ot ti , fa, f 3 , f 4 . 


HvpaTfcft^/. n»m~ 

(a) If |9| < V, 

Shen fj and fa are computed by Eqs. (17) 

fa ■ • ♦ f 1 
f 4 ■ 1 ♦ f a 

(b) If |e| > 4 , 

e 

Sben ccopute f 0 « e 

fa - 1/2 (fo - j 0 ) 

f 4 - 1/2 (f 0 + | o ) 


f x - f a - e 
fa ■ fa - 1 
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HHgtic g**c 


(c) if |e| <3 

Iben fj and fa are computed toy Bqs . (17) 

f a - » - fi 
f 4 ■ 1 • f a 
(a) ife>^ 

Then f* la computed by aaara of Sand polynomials 


fa - 1 - f « 

Also if |8 - aln 0 | < | sin 6 | 

(or approximately If ® < 
then f| la computed by Bqs. (17) 
f 3 - • - f 1 

Otherwise if |0 - sin 8|>| sin 6| 

(or approximately If 9 

then fa Is computed by aeana of Rand polynomials 




*1 ■ 9 - f 3j 

A special problem arises In the computation of the tens Involved 
In Eqs. ($) due to the loss of accuracy In subtracting the nearly equal 
terms Involved* An expression baaed on the binomial expansion removes 
this difficulty* Shis method supplies results more accurate than the 
straightforward computation for terms of the type of: 
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If |B - Bo | Is — U compared to r and Is known mare accurately than can 
be computed by taking the difference between R and Ro« It should be em- 
phasised that this is the case for the Encke term 





since R v , c is computed from 

" *rc + ^ 

and S is small sad known more accurately than can be computed from R^- R^. 

Another ‘"T 1 * will be the sun's and planets' perturbations In 

/ 

Earth reference and the moon's perturbations for a satellite within 3 radii 
of the Earth* For example 

v ‘X* 

R^ is small and known more accurately than can be coaputed from \ B m \ B » 


8pecial Oomnutatlcn of Perturbation Terms 

R - Bo ♦ A 

R Rq R B . B r & 

p • S® • p • to * + v 


or finally: 


-?[ 1 -<!/]* v 
f 1 - Cl * uW a J 


" ro* + 


Bo A . R \ i 

" " V P } i 


i-1 
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where 


u • # <*> * 5 •* 

*» - 'I » jig, •» -^5, 

®» *lx tens are adequate for |u| < .1. for larger Tallies of u 
straightforward computation is adequate. 
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XV. Status of Orbit letexmlnatlcn Progr am To Data. 

Sm orbit determination program, cs it currently exists, le con- 
sidered to be tbs nucleus of vhat is needed for ground based lunar missions. 

In Its final fora the program mil be a complex, general purpose orbit 
d«. arminatioc pragma, requiring the use of a large computer of the LEM 7094 
variety* It Is, however, conceivable that diminutive versions could be evolved 
for special purposes, e.g., saall cooputers at Individual tracking stations, 
and/or spacecraft on-board coaqiuters. 

She current program has been under development for slightly more 
than one year, but la based on experience gained by the personnel Involved 
sines 1956* At that time, the Interplanetary Trajectory Program, based on 
a modified Backs method, was initiated. That program has been widely 
accepted and utilised by M3C and many Apollo bidders and contractors. 

The Backe method Is particularly well suited for computing tra- 
jectories involving the effect of email perturbations acting over long time 
arcs* Sines only ths departure from the local two body orbit is integrated. 

It is possible to compute small perturbations to the full digit accuracy 
machine capacity. Moreover, round-off error due to a large number 
of integration steps, can be effectively eliminated. The Encke method 
divides the desired solution Into two ports; ths first, and larger. Is an 
algebraically determined exact solution of the local Kepler orbit; the 
second, and smaller portion, is the integration of the deviations from 
this Kepler orbit and contains the accumulation of round-off error. By 
suitably rectifying the Kepler orbit whenever the integrated portion be - 
comes too large, it Is possible to control the round-off error. 


a - 



In the process of working with the Interplanetary Trajectory Program, 
the variation of perigee height for the Echo Satellite was computed. 

(Pig. 1). The observed variation. Indicated by the circles. Is compared 
with the computed variations for a period of more than 380 days. This 
single computation carried out for 360 days shows remarkable agreement 
for a numerical program over such an extended number of orbits. This was 
carried out with a fairly crude approximation for the atmospheric model. 

In another example, the inclination as a function of time (Pig. 2) in 
orbit for IMP or Interplanetary Monitoring Probe, computed by two different 
programs is shown. One program developed by Halphen and Improved by 
fir. Risen, and the other the Interplanetary Trajectory Program. This aatell 
la In a highly eccentric orbit and is ctroogly perturbed by the Sun and the 
Mocsi. Halphen 's program *s designed only for the study of the long period 
terms of satellites. Again, the interplanetary Program has shown Its 
ability to retain a reasonable orbit for extended periods of time. However, 
Halphen *a method Is not useful In computing real orbits, because the short 
period terms arc neglected. These results Indicate that the Interplanetary 
Program would be quite capable of producing an orbit determination program. 
The current pregram utilized nary of the features of the earlier program 
and is known as 8PAT-DC for Space Trajectory Differential Corrections. 
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SPAT- DC utilizes the following from the earlier programs: 

A* The modified So eke method— but incorporating a new two body 
solution which eliminates the singularities due to zero ec- 
centricity end zero inclination. 

*• The sixth order backward difference Integration (due to Cowell) . 

C. Planetary Fphnme rides from the Naval Observatory, including 
the Sun, Noon, Jupiter, Mars and Venus. The other planets can 
he readily added. 

The current program is described in Reference 2. It contains a new set 
^Tferantial correction parameters which permits the utilization of obser- 
m ^ ,CB Ofmr long time ares without encountering the danger of a singular 
dtCDu^tlal 'Traction matrix. Reference 3 shows that as sore than qk 
D*»« eter depends on tbs energy, the correction matrix approaches singularity 
**• rapidly. Rjr may of illustration, if one uses the ccmpooents ccf the 
in i t i al position and velocity vectors, Ro and Ro, as orbit differential 
costs ctlao parameters, the solution would not h old up for more th«n jbout 
10 to 12 oxbita* The best that one can hope to do in this regard is to select 
a set of variables in which only one depends on the energy. Such a set Is the 
conventional orbit elements (a, a, 1, «, 0, t p ) . However, these exhibit the 
singularities of ^ero eccentricity# zero inclination, etc. The present pregram 
el i mi n a t es these difficulties, but will require special treatment for the 
Pcxobola and the rectilinear orbit, but with these two exceptions (which are 
•imply added to the existing prcgrasO any orbit that may occur can be handled. 
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To generate the correction Matrix, the program Integrates the 

required Dumber, of separate trajectories simultaneously. It can 

- > 

accomodate up to twenty variables. A conventional least-squares 
technique la then employed to obtain the nev initial Bq and Bo* 
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Variational Psrsnsters 

®> discuss the pirnfri used for variational equation* in the 
Oiffamtlal Correction t tagx m , it la first necessary to review briefly 
tha Algebra of rigid vector rotations. Zhe significant paraaeters in. 
aolvad la rotating vectors are shown la figures 3 and h. Qm U rotations 
are considered, far exa m ple R a and Rj in such a aa inner as to keep t he 
acalar product equal to a constant. 

**> views of tha vector rotation are above. Zhe aide view is 
tfe* plana of the * acton R l and R* and T is the angle throi^h which 
Ra la rotated. 

A coordinate systea is defined such that a a la a unit vector 
B0 * ml a»A In the plane of Rj and R*, e 3 is the unit vector In 

the direction of Rj and eg la the unit vector which defines the right ^ 
ftanded coordinate ays ten. It can he shown that the new vector R* a la ^ 
given hy the equation in Figure 4 . 

9m Initial version of the IASA 8PAT-DC Prcgrm uses seven basic 
variational paraaeters. See figure 5. She flrvt of these, 60, Is a 
mill rotation of the position vector about the velocity vector with 
the resultant position and velocity given by the equations. Shis particular 
rotation has tha effect of changing the axguaent of perigee plus a am . 11 
changes In the I n cl in ation of the orbit and the right ascension of the 
seconding node. 

She second variation, tag Is a saal l rotation of the velocity about 
the position with the resulting aquations shown, and has the effect of 
the Inclination alone, plus saall changes In the azguuant of perigee end the 
right espwnslou of the ascending node. 
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She third variational parameter, t g a is • rotation of the position and 
tha velocity s i mu lt a n eously about the angular momentum. Thus, the entire ellipse 
la rotated about its focus in its own plana. This has the effect of changing 
the argument of perigee. 1 base first three rotations leave the orbit invariant, 
and they only change its orientation in space. These effects are analogous 
to c h a ngin g the I ncl i nation, the argument of perigee and the right ascension 
of the ascending node, without, however, displaying the singularities of 
these latter parameters. 

Sis fourth variation 6o+, is a rotation of the position vector about 
the angular mnmtnlam vector. 

She fifth parameter, 60, is the sole parameter which affects the 
energy in the orbit and is a change in the reciprocal of the seml-Mjor axis. 

Ihie is done by also changing R and R in such a way as to keep the eccentricity 
equal to a constant. 

She sixth variation, 60+ is a change in the Mgnltude of the position 
vector. As a result, there must be a corresponding change in the 
of the velocity vector to keep the energy invariant plus a rotation of the 
position vector about the angular momentum vector to keep 60+ invariant. 

She seventh and final variable, 607 is a snail change in the drag co- 
efficient. She parti ale are then obtained by integrating seven different 
trajectories each corresponding to a particular tp. and the coordinates of the 
computed position subtracted fros> the coordinates of the reference trajectory 
and divided by the variation in a secant approximation to the derivatives. 

Qiese then are the basic constituents of the program. 

SO reiterate, the main advantages are that the parameters are chosen 
so that only one depends on the energy end ell singularities have been 
eliminated. 
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VI. OtLfflereotlal Correction Methods 


Since the orbit position and velocity ie not directly observable, it is 
necessary to infer these variables from & sequence of observations which are 
functions of the trajectory. Zh the conventional methods, a linear relationship 
Is assumed between the deviations In the observations and the corresponding 
deviations In the orbit variables. Thus, an error In the orbit position will 
correspond to some predictable error In the observation. A large number of 
observations are made, overdete raining the linear system of equations, a 
least square technique is used to obtain the best value of the orbit errors 
to fit the known observation errors* Since the equations of motion axe essentially 
n o nline ar, this linear region becomes more and mare constrictive about the 
no m in al trajectory the longer the time period over which the prediction is mede. 
Thus, the least square technique often produces a result, fittlr^ data over a 
long time arc, which is outside the linear range. This produces problems in 
convergence and consumes machine time. Reducing the number of observations 
to a shorter time arc helps avoid this difficulty. The weighted least squares 
is often used in this manner* However, a large number of observations Is 
always needed in order to properly evaluate the effect of the random instrument 
errors. 

A. Minimum Variance 

The Schaddt-Kalman 'deferences 4, 5) minimum variance technique 
avoids the difficulties mentioned above. This procedure permits a complete 
optima estimate of the orbit variables and the observation errors from each 
si n gl e observation. It has been adequately shown that the two methods converge 
to the same answer eventually for the same total set of data (Reference 6). 

Moreover, the variance technique always converges more rapidly since it requires 
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less data at each stage than the least square procedure. At best, the least 
square technique can, he said to he as good as the minimum variance* 

It Is vorth briefly stating some of the Important formulae for the 
mi nimu m variance technique. The nominal trajectory may he obtained by 
integrating a system of differential equations with certain unknown parameters, 

x i “ f i * x j * a • • (20) 


Since these parameters are in general not observable. It Is 
necessary to make observations whose predicted value la given as a function 
of the trajectory coordinates 


y - y (x , t) 



( 21 ) 


The deviation of a vehicle from its nominal Instantaneous position, 
due to a change In the orbit parameters, may be obtained from the linear term 
of the Taylor series 


X d x, (t) 

4 *i (t) * ^ 5T ; ( 6 ) tx i (0) • (22) 

The matrix of partial derivatives (*) Is called the state transition 
matrix and relates the expected deviation from the nominal trajectory at some time 
(t) vith the deviation at seme previous time (t Q ) 

A \ (t) - * (t , t Q ) d x t (t Q ) . (23) 
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An estimate of the predicted deviation In the observotlor 
corresponding to the deviation from the nominal trajectory, Including 


obse r v at ion error s , is given by the linear eqoation 

A y (t) ■ Z j x^^ ft) A * ± (t) + « • (24) 

2be covariance matrix of the squares of the expected deviations 
of the octlt variables may be derived from Bq. (23) 

P (t) - t (t , t Q ) P (t Q ) t* (t , t Q ) . (25) 

Uhls equation prellcta the change In expected error In the orbit 
as tine Increases* The covariance matrix of the observations as derived from 
Bq* (24) la given by 


X (t) - M P (t) M* + 1 a 

" (t) ■ (H-$) . 


( 26 ) 


It is required to find the optimum estimate of the relationship 
between the orbit errors and the Observation errors In a form inverse to 3q* (24)* 
The optimum estimate for the solution of Bq* (24) from a single 
observation Is given 


A x (t) - K (t) A y (t) 


(29) 
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Estimating A y from it* expected value in Eq. (24) 

A x (t) ~ K (t) C M A x ♦ « ^ • (30) 

The equivalent covariance equality for Eq* (JO) is given by 


PM* - K (MPM»+? ) . 


(51) 


In the above, it is assumed that the observation errors are 
uncorrclctcd vith the orbit position errors. Solving the Eq. (31) for the 
optimum estimator K (t) end substituting into Eq. (29) 

A x (t) - P (t) M* (t) Y~ x (t) A y (t) (32) 

y 

This solution moy he used to correct the orbit positions eh the 
end of each observation* Furthermore 9 the covi-riance matrix of the orbit 
errors may be reduced to a smaller uncertainty following each observation 

P + (t) • P* (t) - P~ (t) M* (t) Y* x (t) M (t) P" (t) . (JJ) 
where P” (t) is the value of P (t) obtained from Eq. (23) 

Using the best prediction of the orbit position, a new trajectory 
may be integrated following each Observation. In addition, a new stats 
transition matrix and covariance matrix way be computed to predict the next 
observation. In this manner, the complete orbit information is accumulated 
about a nominal trajectory of higher and higher accuracy. 
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B. Weighted Least 8quares 

The aethod of least oqusrcs any also be used to obtain a solution 
of the deviations in the state variables corresponding to the actual 
de vi ations in the observations* Equation (2b) any be re written to relate 
the deviations in the observations to the deviations In the state variables 
at the initial tine. 


A y (t) - N (t) t (t , t Q ) A x (t Q ) + c • (54) 


Shis aquation nay be written for an entire sequence of observations In order 
to font an ovexdetendned aysten of linear equations for the unknown A x (t Q ). 

The solution of Eq* (3b) by the method of least squares Is given by 


A x ± (t Q ) - (A* , A)" 1 A* A y 


there A ■ N (t) f (t , t Q ) and where A y is the actual deviation in 
the observations* 

Ih this solution no a priori knowledge of the expected deviations In the 
observation Is required or assumed known* 

The value of the Initial parameter Obtained from Eq* (35) nay be 
used to derive a new nominal trajectory and a new set of residuals, to obtain an 
improved estimate of the Initial state variables, A x (t Q )* The process may 
be repeated until no further improvement is obtainable from the observation 
data* 

In order to speed up the rate of convergence, it Is necessary to 
I n clu d e estimates of the expected deviations In the observation through the 
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use of covariance matrices and the method of weighted least squares. The 
solution of the Initial state variables by the method of weighted least 
squares is given by 


A \ (t©) - [ A* V A + ]' 1 A*W4y (36) 

where V - (l*)* 1 


P i-l “ Cov [ A x 1-x (t o ), A x 1-x <t Q ) 3 


The essential difference between the method of minimum var iance , 
least squares and weighted least squares is contained in a comparison of 
equations (32), (35), *nd (36). 

The method of least squares and weighted least squares both relate 
the estimate of the initial parameters to an entire sequence of observational 
residuals spread over an extended time arc. In contrast, the ernthod of 
variance relates the present estimate of the state variable deviation to the prese: 
actual deviations in the observations. The linear assumptions required for the 
updating theory are violated to a much less degree in the method of minimum 
variance than in the method of weighted least squares. 


- 32 - 



VII. InOntiQB of tbs Partial Derivatives 

R Is seen th a t both the least square technique and the minimum 
iraslaBea technique require the computation of the partial derivatives of 
the observations with respect to the variational parameters. As one 
*** tj ro * c b these sear be obtelned by integrating additional trajectories and 
foodag the differences using the secsnt method. However, the program 
xecCBMeaded in this report obtains the matrices of partial derivatives 
analytloallor In terms of the Instantaneous two body orbital coordinates. 
Xhus # the complete orbit prediction and partial derivative matrices may 
be obtained in ess entiall y the same computing time as that of the nanl nal 
trajectory. 

2he partial derivatives are obtained from the product of the M 
discussed In Section VI, and a state transition matrix. The 
method of obtaining the state translation matrix Is based on a generalisation 
of an tncke method applied to linear prediction theory. It Is assumed that 
*be sqpiatlona of motion may be decomposed into two factors 

i 

x - g (x , t) ♦ b (x , t) ( 58 ) 

where 


I S | > > | h | . (39) 

It Is further assumed that a closed fom solution of the differential 
equations is known for the case vben b ■ 0 , 

• “ g (s , t) . (bo) 

furthermore, the state transition aatrlx for the app r oximating solution 
la known In dosed form 


As (t) - * (t , t o ) &s (t 0 ) 


( 41 ) 
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Let the deviation between the state variable and Its approximation 
be given by 

P (t) . x (t) - • (t) (kl) 

®je perturbation eolations of notion may nov be written In the 


e 

P « g (x, t) - g (s, t) 4 h (x, t) . ( 1 * 2 ) 

In order to guarantee that the deviation, p, la never permitted to grow 
too large, the process of rectification Is Introduced. Whenever a pre- 
determined value of p is exceeded, the Integration is terminated at timt 
* r * * *•* of initial conditions are introduced, setting p (t ) equal 
soro* Integration proceeds again about this aw appro xl rmts 

solution. y 

8lnoe the deviation between x and s is never permitted to exceed 
the given value, the partial derivatives of the state variables from 


their nomin a l value may also be limited, 
an approximate state transition matrix 
< (t, t o ) r , ( t , t<> ) 

for t s t. 

r 


Ihus It is possible to write 


( 1 * 3 ) 


Moreover, the approximate state transition amtrix Is known In closed 
form frcm the knowledge of toe known solution of the trajectory 

for Eq. (4l). Following each rectification, it is necessary to relate 
the state transition matrix at time t to the time. Bila may be 

accomplished by miltlplylmg the approximate state transition matrix 
for times within each rectification Interval by Its value at toe last 
rectification. 
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# K * 0 ) 2* f (t, t 0 ) k f (t, t ) 4 (t r ^, V (44) 

Xft this manner, s closed fora expression for the state transition 
Matrix asy tea Obtained without Integrating large quantities of differential 
stations* Moreover, the error in the computation nay be limited by 
the proper use of the rectification technique. 

As has been stated, the use of the conventional state variables, 
nanely initi a l position and velocity, leads to a transition matrix which is 
poorly c ond i ti oned and Uhlch contains rapidly varying functions of time 
for its clement a. It the procedure recommended here, the first six 
variational parameters described ?- Section V have only one secular term, 
namely that due to a* 

lb addition they have the characteristic that they completely 
determine the orbit Independently vf the orbit orientation or shape and 
never break down. 

A co m plet e modification of the Schmidt-4(alman equations in terms 
these new parameters has been derived and is now available for 
incorporation In the prediction scheme. 

lbs n od i fle d Schnldt-Kslwan equations are listed below. Deviation 
of the orbit variables Interns of the net. parameters are given by 

A* (t) = f j & (t) S (t, t) &.(*) 

' * ' (45) 

0 (t) = M S (t) 3 H (t) to (t) . 

®xe nev transition matrix It given by 

A® (t) = C ) Aa(0) = o(t, t Q ) rjx (0) . (46) 


( 
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The observation errors in terns of the new parameters are given by 
Sf \t) B l(t) faX (t) + ( 

( 47 ) 

vhere H(t) ■ M(t) 8 (t, t) . 

The corresponding covariance matrices are given by 

Q (t) , n Q (0) cf 

P(t).8Q (t) 8* (48) 

y (t) xr h q (t) ** + j* 

The inverse relationship between the orbit parameters and the 
observational errors is given by 

da(t) * L (t) sf (t) . (49) 

The optleue estimator L (t) is given by 

L (t) = Q H* y" 1 . (50) 

The corrected estimator after each observation is given by 

Q + = <f - Q’ J* H Q* . (51) 

Using these modified equations, it is now possible to 
use the Bo h m l d .t - K a l m e n scheme far bath short and long t era predictions. 
Moreover, the computing time need not exceed that necessary for a 
single ncmlnal trajectory. 

It is sometimes necessary to obtcin an estimate of the 
instrument errors from the accumulated experience of many observations. 
Let the observed error between the predicted and actual observations be 
given by (y). Then, an estimate of the instrument errors, after many 
observations, nay be obtained by 

•* (t) » (I - I L) (y) (I - H L)* . (52) 
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*» lUuttnt* tba (Illicit; «4 brant; of tb*M 

taUani * “W"" of jmm dralratlra. « 






«bm 


Qm m f*ul { ( -Uj- *.*,) } 

' 1 r*(t) v h* J 

and 

■*^v u {•ffo [ f * l,) +_ fto J *^r * 

♦2s u -‘’w^ rt] ) 

9mm mtw to Indicate how simply all the required aatrices 
■ay be computed as functions of the noainal trajectory. The quantities 
«* and idilch appear in the aatrlx Q are defined as 
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VIII. types of Observations 

One of the advantages of the Schmidt-Kaiman scheme is the ability 
to determine the optimum observation to be made at each orbit position 
in order to have the greatest decrease in uncertainty. The method permits 
a predicted estimate of the decrease in uncertainty for each observation 
prior to having made it. In this manner, a choice may be made to obtain 
more rapid convergence to the proper solution. It Is possible to use 
range and/or raigte rate at each position of the orbit in order to obtain 
the aavlmnm information to be gained by using each or both. In addition, 
observations made from the vehicle during aid-course, from accelerometers, 
optical measurements, etc* may be used in an interspersed manner to optimum 
advantage. It has been known for some time that the simplicity in design- 
ing, constructing and using c v Doppler (Reference 7) range and range rate 
Instrumentation greatly favors these types of observations. 

The 8cfamldt-Kalman technique, by procc-sslrg the data within the 
linear range of the prediction theory, and through the use of the recommend- 
ed orbit correction parameters, makes it possible to rely on these instru- 
mentation techniques over large portions of the tracking complex. 

A. Transformation Equations for Tracking Data. 

It is sometimes desirable to transform the computed data into a 
topocentric coordinate system. The computed data is obtained in a cartesian 
coordinate system (x,y,z) defined with the origin at the center of the Earth, 
the x-axls in the direction of the vernal equinox and the x-y plane in the 
equatorial plane of the Earth. The series of coordinate transformations 
required Is as follows: 


A rotation about the z axis through angle G' so thac the x ' -z ' 
plane is in the Meridian plane of the stations 
Then: 

x* • x cos 0 ' + y sin 0 ' 

y* ■ -x sin f + y cos 0 ' (5 1 *) 

s' ■ X 

where 


x y s' are the rotated coordinates and 0 ' is the right 
ascension of the Meridian at the tine of observation, e* is computed 
by adding the longitude of the station to the right ascension of Green- 
wich at time of the observation. The right ascension of Greenwich for 
the beginning of the year in question is obtained by Newcomb's formula, 
updated to a i960 origin. The instantaneous hour angle is then obtained 
by a linear formula, separatirg daily and hourly rates. 

A translation of the x', y z * axes from the center of the Barth 
to the station in question 

x" • x' -(C * h)'cos ® 

y" » y' (55) 

z" » z' -(S 4 h) sin <p 

where 9 is the geodetic latitude and h the height above sea level cf 
the station in question. C and S are parameters which account for the 


polar ohlateness of the Earth and. are given by 

C • (1 - a* sin* 

8 - (1 - a®) C (56) 

A rotation about the y* axis through an angle 90° -<p so that 
the a* axis la in the senlth 

x'" a x" ain <p - t" cos <p 

y- - y* ( 57 ) 

s'" ■ x" cos ® + z" sin <p 

Shis series of transformations will transfora the geocentric 
coordinates of a point into the topocentrlc system. Transformation of 
topocentrlc to geocentric coordinates is achieved by the following set: 

x" « x'" sin ® * s'" cos a 

y" ■ y (56) 

1 " * -x cos 9 + s'" sin © 

x' » x" + (C + h) cos 0 

y* - y" (59) 

t' ■ z" + (S + h) sin <p 
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x ■ x* cos 0' - y' sin 6’ 

y - x' »in 0' + y* cos O' (60) 

X * s' 

It is dear, that application of these transformations requires 
the complete knowledge of the topoceniric displacement vector. The 
following paragraphs describe the treatment of various sets of Incomplete 
data that nay be available. 

B. Formulae for observed variables and their Partial Derivatives . 
The progr am will accept the following types of observational data, 
singly or in combination: 

1. Range 

2. Range rate 

3. Right ascension and declination 

4. Azimuth and elevation and Minitrack observations 
In order to generate the differential corrections described in 

Section VI, it is necessary to compute residuals which consist of the 
difference between computed values of the observables and the observation 
data. In addition, it is necessary to compute partial derivatives of the 
observables with respect to the orbit parameters. The range, range rate, 
right ascension and declination can be expressed directly ir terms of the 
geocentric state variables and the required partial derivatives may be 
obtained as follows: 

Range: The computed value of the range is given by 

p •*[* ( x - x B ) 2 + ( y - y 6 ) 2 + (t - * s ) 2 1 ^ (6D 
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■®d the Matrix of partial derivatives for the range is given by 


#t) - [ 



HI* 

p 


* - x 


0 , 0 , 0 


1 


(62) 


Bate: She computed range rate is given by 


P " “ f U - x g ) (* + «> e y B ) + (y - y g ) (f - ® e * B ) + U - x g ) 


1 


(63) 


She Matrix of partial derivatives of the range rate with respect to the 
state variables Is given by 


*tl.r * * "• — _ • $ * ®e ** P (y-ys ) £ p (*-*a) 

* o p T • Z — > — 


2 P 2 

P P 


(64) 


x - *s * - * 


1 1 


Ascension and Declination: The expressions for the right ascension 
and the declination nay be written as 


sin D 



tan BA 



(64-a) 
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itrix of partial derivatives Is given for D by 


x - x y - y (*-*,)*“ P 2 

* T’c'STd ! — ~ * “ — ' xzzrt ' 0, 0, 0 1 


p ' pU-* g 5 


and for R A by 


1 r y ‘ y a 1 1 

*TG k ! * TTT , o, o, o, o I 


sec 4 RA 


(x - a ^ 2 


X - x_ 


Azinuth and elevation and the Mini track observations are nost 
conveniently expressed In a topocentrlc, local horizon coordinate systen 
and to treat then It Is useful to Introduce the follovlxg relation between 
the topocentrlc and geocentric coordinates: 


I sin - cos S 1 


- sin V 


sin s sin ?’ - cos 


cos S’ 


9 I [*‘*1 


i J j cos •• cos 9* cos c *ln 5' sin © I I z-z 8 


This relation is used in developing the required partial derivatives 
for these angular observations. 

Azinuth and Elevations: The expressions for azinuth and elevation are 


- 1 y'" 

A » tan “th 


E » tan 


<o 2 -*'" 2 ) 


i» 2 \i 
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Am corresponding attzieu of partial derivatives are 


for A 


tft) 


f-**" *ln e'-y'" cos B' ala 9, 

ft -«••* 


cos e’-y'" aln 0* *la 9, y*" cos ^ o, 0, o"| 


( 67 ) 


*t) 


( •• cos 9 - ~ (x-x^), 

/,. ’ 

•la 0* cos^p - (y-y ), 

P a 

•la 9 - 5 -— 0, 0, o "| 


( 68 ) 


Mini track Observations: 9 m Mini track syntax direction cosines are expressed 
tarae of the topocantrlc coordinates as 


/ - - 


r»* 


111 

P 
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n ■ 

P 



tto oorraapoodlag aatrleaa of partial darlvatlwa or* for 4 


. , 1 T x , *(x-x ) *'"(y-y ) 

*(*) ■ Jl-ntof ooi ♦ , -ala f ala •• A £ 

' o* »* 


*•*(»-»,) T 

« 0 » f ♦ p 5 — , 0, 0, 0 I 

p 


for a 



M(t) 



•la •* 


y* *(»■«,) 

“7 


ooa I* 


y*"(y»y,) 


y # "(a-a g ) 

~ 


* Of 0, 


0 


] 
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and for a 


Mtt) - i [ooa ? coo *• - 


• •*(x-x ) 


, coa ^ ala 6' - 


*'"(y-y.) 


•la 9 


i»"i 
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IX. Current Work 

9m discussion Milch fallow outlines sane of the work to he dene 
over the next >4 year period. 

X. fits current p ro gre s s uses the conventional least squares correction 
techalqM. Xt varies,, and has been tine tested, but there are better and 
■ore Modern techniques, such as the Kalaan-Sdmddt which has shown promise 
of being far superior. The major reason for this is that both techniques 
are using linear theory to correct non-linear effects. However, the conventional 
least squares requires an over-determined solution which in real 

lift means male* eg to Obtain a large nusber of observations during which 
tine the non-linear effects build up. Conversely, the Kalman-Schmldt uses 
each observation as it occurs and thus is ouch closer to the linear range 
on which both theories are based/Another is the one currently In use at 
JEL, which briefly is a weighted least squares technique. Therefore, one 
of the first problems is to develop the Kalaan-Schmidt, and the weighted 
least squares techniques for the present program, to thoroughly test them, 
and compare the techniques. 

2. During the checkout of the present program, the problem inherent 
in many orbit determination programs to date has also occurred. Namely, 
the problem of Ill-conditioned data. Furthermore, the problem varies with 
the type of orbit involved, e.g., for highly eccentric orbits angle data 
alone does not determine the energy very veil, but for nearly circular orbits 
they do not determine the orbit orientation in space very well. For range 
and range-rate data, the converse is true. . As a result, the problem of 
increasing the rate of c o nver gen ce or pursuing several other possible avenues 
to overcame the ill-conditioned data problem requires extensive research 
and development. 
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3 . Many time details have not been ccnpletely accomplished, 

e.g., atmospheric Models, ionospheric refraction, nutation, lunar llbratlan, 
choice of reference equinox, selection of optimum time standards, adequate 
conversion from ephemerls time to the appropriate time standard to name Just 
a few. 

4 . It is recognized that Integrating the number of trajectories 
corresponding to the number of variables is a poor and time consigning way 
of obtaining the correction matrix of partial derivatives. It is felt that 
here again the Encke techniques will offer great advantages over other orbit 
determination techniques, if one can make use of the analytic two body orbit 
to determine the derivatives. If this works, we will be bade to integrating 
only one trajectory. Ihere are, however, problems — and research and 
develop nt In this area is required. 

5 * Before the program can approach anything like a real time program, 
the entire progr am must be examined w.r.t. time consumed In each operation, 
each sequence of operations and each subroutine within the confines of the 
required or desired accuracy* 2h addition, the program should automatically 
select its own optimum integration interval — but here, again, is a research 
and development study since there is a trade-off on how long it takes to 
change integration Intervals, and how much is gained by affecting the change. 
These selections may also be orbit type sensitive. 

6 . The final program also make use of spacecraft on-board 

observations, optionally, alone or in conjunction with ground observations. 
Therefore, the compatibility of the two systems should be coordinated. In 
addition, the output fraa this program should bo useable, when desired, to 
up-date the on-board information. 
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7 . Ibex* tre • myriad of detailed arrangements to staaiardize in- 
obsarvaticnal data to keyed formats ae veil as outputs* Suitable 

ecnieeilflM to 1 he teat network of station a etc. 

8 • 9m program s h ould be adaptable to the rendezvous mission and It 
■V ^ d esirab l e not to oamslt to a pre •designed trmnslunar trajectory, but 
rather to dateralne the appropriate one after the rendezvous has been afft'rfced. 

Xu fhetf la light erf the llalt l e aa amber of possible lunar trajectories. It 
W desirable to select the traaslunar trajectory, even for e single booster 
launch, d> 1 3 « the spacecraft is in the parking orbit. 

9 . lhere are positive indications that it nay be possible, through the 
Ka lm a n-flc famldt technics to design a program capable of differentially correct- 
® rai 9 m ground up, i.e., during boost. If this were so, one of the 
largest dlffl oil ties of all orbit determination programs would be eliminated, 
i.e., the error ia beginning vith a preconceived (nominal) R and V. This 
error is largely responsible for such tracking errors as wrong lobe Identification 
from tracking stations etc. 

10. Finally, and certainly one of the nost important, is the problem 
of bui l d in g in the decision making into the program wherever possible. Bi<« 
F«*lem is somewhat Interdependent with number 9 , since the type of decision 
diich nast be made depends on how large an error must be accounted for. Pii* 
is particularly true for linear theories such as these because any malfunction 
can cause the solutions to be cutslde the linear range. The present program 
rejects all observations which are outside the 3 c deviation for each Iteration, 
but this only scratches the surface of »»rvn 4 ng the problem of bad data. JPL 
baa a very extensive data editing prograa. Related to this is a study of 
maximizing Information returned from various types of tracking data, e.g., 

vange data, ibich may very well dete inine the energy of hi g hl y eccentric 
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orbit*. ** ■y Possible to neglect other variables and get a first order 
approximation of drag correction, similarly refraction corrections etc. 

11 . In all this, &ev concepts, variables and programs are being intro- 
duced* Zt s hcwld go without saying that a continuous and exhaustive shake- 
down of all cf these under all conceivable conditions ahmii* be a ecus tent 
part of this total effort, e.g., determining the characteristics of the 
differential parameters under all types of ozblts, at various positions in 
each of these orbits, with ideal observational data, noisy generated 
observational data and wherever possible, actual satellite or apace probe 
data* 


It also goes without saying that all this work be properly 

documented with progress reports, interim and final reports. Puoerous 
prcblcsu will also arise during the course of this work for which solutions 
■ist be found to the mutual satisfaction of all. 



X. Ba fraction Correction 


A. Refraction Correction for Pulse Radar 

B the event that the observation signal Is an actual measurement 
of a time delay la. transmitting a radar pulse, the assumption is usually amdc 
that the rango la obtained bgr simply multiplying the transit time by the 
vacuus speed of light, c, as follows (see figure 6 ) 


«et 


. f 


r v c 

- d s 
r o v 


( 71 ) 


lh term of the In d ex of refraction, n, the apparent, range is 
» j r n ds . 


( 72 ) 


She actual path traversed by the signal is given by the rule that the transit 
tlae is a a lrrtiaim . From the calculus of variations, the solution of the path 
Is characterized by the fact that 


k • constant m n(r) r cosy * n(o) r Q cos . (73) 


Be f erring to fig. 7 the relation between the arc lex^th, ds, and the path 
angle, Y, is given by 

n r dr 

" Jn 2 r 2 -* 2 * (7^) 

Equation (72 ) may be written 


p. - 1 ry 

* J r o Vn 2 ^-k 2 


( 75 ) 
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Eq, (75) my be Integrated in dlcsed fom provided n(r) is a constant 


P 


i+l 






k 2 - v n 2 xf - k 2 


( 76 ) 


Assucii^ that n(r) la piecewise constant over several layers in the 

atmosphere, It Is possible to obtain the apparent range by summing the successive 

solutions of Eq* (76) until the cosiputed geocentric distance r y is obtained. 

Since the constant, k, is given in terms of the initial sight angle, 

k ■ n(o) r(o) cos X , the apparent range is seen to be a function of the 

o 

geocentric distance, r y , and the grounl station sight angle, X Q . 


P. ■ Z (p. - P* ) 
a i-1 * 1-1 


( 77 ) 




E. Refraction Correction for Range or Range Rate of CV Doppler Radar 


Data. 


When the signal observations for range or range rate are obtained 
from a CV carrier measuring either phase shift or frequency shift, it is 
necessary to correct for the additional change in frequency or phase due to the 
effect of the refraction on the phase velocity of the signal. A first order 
estimate of this correction may be obtained by using the concept of the average 
index of refraction over the signal path. The Doppler shift Is given by 


A f 
* R 


( 78 ) 


where A f is the difference in frequency between the emitted and returned 
signals, and n is the value of the average index of refraction defined as follows 




n » 



(79) 


(79) «q r Integrated wing tbs aeeusptian that the index of refraction 

la placentae oonetaat 
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£ “i (h l+l * 

v-*; 



( 80 ) 


®» al titu d e of the vehicle my be obtained fron the aaamed comuted 

poaltlon of the vehicle abore the geold* To apply the refraction correction 
for G V p haee range and range rate observation* It la only neceseaxy to replace 
the mine of the speed of light la vacuus by c/n In the noonaal equations used 
to obtain range and range rate eetlnates frcn the signal* 
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XI. Rec c — c ndu t loos 


It la recoansnded that a general purpose orbit prediction and 
tracking program he generated incorporating the features outlined in 
this report. 

It la further ra n n— ende d that specific attention be paid to 
the requirements for the Apollo space mission so that this prograa any 
be used la the Apollo tracking oonplex, or as a back-up for whatever 
system is finally decided upon. 

It is recoameaded that a detailed comparison between the methods 
suggested herein end the least square procedure, using other parameters, 
as well as the weighted least square procedure recommended by J.P.L. be 
extensively explored so that the best features jof each may be Incorporated 
lx> the final acceptable scheme. 
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RIGID VECTOR ROTATIONS (1) 
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RIGID VECTOR ROTATIONS 
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M«ure 5 (2) 



Th* final R and R fear all alx £t'a la than 



Sba last parwatar change does not affect tha Initial condition# but 
entara Into the acceleration confutation. 
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Figura 6— Typical ray path 





